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Abstract. The three-dimensional baroclinic quasigeostrophic flow model has 
been widely used to study basic mechanisms in oceanic flows and climate dy- 
namics. In this paper, wc consider this flow model under random wind forcing 
and time-periodic fluctuations on fluid boundary (the interface between the 
oceans and the atmosphere). The time-periodic fluctuations are due to pe- 
riodic rotation of the earth and thus periodic exposure of the earth to the 
solar radiation. After establishing the well-posedness of the baroclinic quasi- 
geostrophic flow model in the state space, we demonstrate the existence of the 
random attractors, again in the state space. We also discuss the relevance of 
our result to climate modeling. 
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1. Introduction 

We will study the quasigeostrophic equation (QGE) for the ocean in a three di- 
mensional domain. This geophysical flow model has been formally derived as an 
approximation of the rotating three-dimensional primitive equations [ p^ . Bour- 
geois and Beale Embid and Majda [|l^, and Desjardins and Grenier Q have 
recently shown that the three-dimensional quasigeostrophic equation is a valid ap- 
proximation of the primitive equations in the limit of zero Rossby number. Holm 
|l4| has established the Hamiltonian formulation for the inviscid quasigeostrophic 
equation. In particular, we will consider the following version | pl| , |^] of the equation 
in terms of the stream function ipi^i Ui t)- 

( P 

with AV' = -^xx + V'yy + ( jy2°^-) ^^ 

Here a;, y, z are Cartesian coordinates in zonal (east), meridional (north), vertical 
directions, respectively; /o -|- /3y (with /o,/3 constants) is the Coriolis parameter; 
N{z) > is the Brunt- Vaisala stratification frequency, and v > is viscosity. More- 
over, J(/, g) — fxgy — fyQx is the Jacobi operator and potential vorticity is defined 
as A-f/) -t- /o + fiy- Note that A-0 can be regarded as a modified Laplacian operator 
where the coefficient in the vertical z direction is adjusted due to the density strat- 
ification, and the coefficients in x, y directions are constants due to the horizontal 
density homogenity in the 3D quasigeostrophic flow model formulation. Bennett 
and Kloeden have also used a similar modified Laplacian viscous term in a more 
complicated 3D quasigeostrophic flow model involving thermodynamics as well as 
hydrodynamics. 

Our aim is to study the potential verticity evolution in an ocean under the influ- 
ence of the atmosphere. Following an idea of Hasselmann one can divide the 
geophysical or meteorological flow into two parts. These two parts are the slowly 
changing climate part and rapidly changing weather part. The weather part can be 
modeled by a stochastic process such as white noise, see Hasselmann Arnold 
|2| and Saltzman [^. Oceanic flows are affected (on the ocean surface) by these 
short time influences due to weather variations which are usually called wind forc- 
ing. Moreover, oceanic flows are also affected by climatic variations due to periodic 
rotation of the earth and thus periodic exposure of the earth to the solar radiation; 
see [|0|, Chapter 6 and [|6|, Chapter 11. 

Since the exchange between the atmosphere and an ocean takes place at the surface 
of the ocean, we will consider the above quasigeostrophic partial differential equa- 
tion with white noise Neumann boundary condition on the top surface of the ocean 
|T9[ . We will also add some time-periodic boundary condition on the ocean surface 
due to periodic rotation of the earth and the solar radiation. Since there is no 
inffuence of the weather at the bottom of the ocean we will have the homogeneous 
boundary condition there. The boundary conditions in horizontal directions are as- 
sumed to be periodic for mathematical convenience as in other recent works ^ . 
However, since we have a forth order differential operator at the right hand side of 
the quasigeostrophic equation, we need a second group of boundary conditions, say, 
for Kij}, and these conditions are the ones used in |^ to be specified below. 
For simplicity we will suppose that this second group of boundary conditions are 
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deterministic. But generalizations are possible. 

Our aim is to find structures in the dynamics of the QGE. As we have time de- 
pendent random and time-periodic boundary conditions, we will obtain a nonau- 
tonomous dynamical system with random influences. We will show how to find 
attractors for such a dynamical system. The existence and interpretation of cli- 
matic attractors have been controversial and have caused a lot of debate Q. A 
low dimensional climatic attractor was regarded as an indication that the main 
feature of long-time climatic evolution may be viewed as the manifestation of a 
deterministic dynamics. Our result is about random attractors, and thus the long 
time regimes that such attractors may represent still carry the stochastic informa- 
tion of the geophysical flow system. Stochastically forced QGE has been used to 
investigate various mechanisims in geophysical flows |2^, ^ . 

The QGE can be transformed into an evolution equation with standard bound- 
ary conditions. For this transformation we need an Ornstein-Uhlenbeck process 
fulfilling our dynamical random or time-periodic boundary conditions. This trans- 
formation will be introduced in Section 2. 

In Section 3 we investigate the coefficients of the transformed evolution equation, 
and further obtain a global existence and uniqueness result and some regularity 
result. 

In Section 4, we study the random dynamics of the transformed QGE. Based on 
the uniqueness result above, the transformed evolution equation generates a nonau- 
tonomous dynamical system. In addition, if we restrict this system to discrete time 
step of the period of the periodic rotation of the earth, we obtain a random dynam- 
ical system. This random dynamical system has a random attractor. This result 
can be extended to the dynamical system on the real time axis. 
The Section 5 contains the proofs. Finally, we summarize our results in Section 6. 

2. The 3D quasigeostrophic equation 

Let O = (0, 2tt)^ be the cube which is a model for a piece of the ocean. For 
X, y, z G (0,27r) and smooth functions u{x,y,z), v{x,y,z) we define the Jacobi 
operator 

J(m, v) = UxVy - UyV^. 

In addition, the differential operator A is defined by 

A:=d.,^+dyy + iF{z)d,)^ 

for /o £ M \ 0, F{z) := j^ij^ and N{z) is defined to be a positive C°°-smooth 
function N{z) > on [0, 27r]. Let i^, (3 be positive constants. In the following we 
investigate the 3D QGE flow model ||, H^]: 

(1) A^Pt + J{i^, Aij) + 134,^ = i/AAV'. 

We impose the following boundary conditions for this equation. Let O-,-,0j 0.^.^2tt, • • • 
be the faces of the cube O then we assume periodic boundary conditions 'm x, y 
directions 
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and similar for the faces O..o,-, 0..27r,- together with derivative in y direction. With 
respect to the z direction, ip fulfills Neumann boundary conditions on top and 
bottom boundaries; see p. 367. Namely, we impose that 

dtp 

— — — — Q on bottom boundary O. . o 
an 



and 



dtj} 

-77- — 'ipz =W on top boundary O. . 2-w 
on 



o 



where is a temporal white noise and n denotes the unit outer normal vector. 
Moreover, Aip is supposed to be periodic in x, y directions: 

and similar for the faces O.,o.-j 0.^2tt,- with derivative in y direction. For the other 
faces on the top and the bottom of the ocean, we impose homogeneous Neumann 
boundary conditions as in [^l], ^, 0] 

AV'.|a,.,o =0, AV'.lo,,., =0 

In addition, we impose that 

/ ipdO = 0, 
Jo 

and 

AtpdO = 0. 

We also assume an appropriate initial condition 

ip{x,y,z,0) = tpo{x,y,z). 

Later on we will see that we can find an Ornstein-Uhlenbeck process rj satisfying 
the linear differential equation 

(2) Vt = i^A??, 

where 77 fulfills the same boundary conditions as ip in some sense. 
We now transform formally the above equation to another parabolic differential 
equation with random coefficients but with standard boundary conditions. To do 
this we have to introduce the solution operator G of the following elliptic problem: 

(3) AC = / 

on O with periodic boundary conditions in x,y directions and homogeneous Neu- 
mann boundary conditions in z direction. For existence and regularity properties 
of this operator see below. G can be considered as a linear continuous operator. 
Here / is an element in a particular Hilbert space. 
For the difference of ip and rj we obtain for 

{Alp - A7])t + J{ip - i],Aip - Ai]) + J{i],Atp ~ Aij) + J{ip - 7],A7]) + J{7],A'q) 
+ P{ip~r])^+pij^:^i^AA{i;-r]) 
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where 77 is assumed to be sufficiently regular. The difference ip — r] satisfies the same 
boundary condition as ^ in (^) . Hence by the formal substitution v — At/j — Ary we 
can write 

vt + J{G{v),v) + J{ri, v) + J{G{v), A77) + J{ri, Kif) 
+ l3G{v)x + I3r^x = vAv. 

The problem (|2|) satisfies only standard boundary conditions with respect to 77 but 
not with respect to A?]. To obtain an equation with standard boundary conditions 
we need the transformation u — v + Arj assuming formally that 77 is sufficiently 
regular. Then u fulfills periodic boundary conditions in x, y directions, and satisfies 
the following random evolution: 

(4) ut + J(G(w), u) + J{t] - G(Ary), u) + pG{u), - P{G{Ari)), - 77,) = vKu. 
This equation will be treated in the rest of the paper. 

Remark 2.1. To obtain a solution of the original equation all reverse transfor- 
mations make sense if 77 is contained in a Sobolev space of second order. Hence 
we can find a solution of this equation contained in this Sobolev space. Thus, in 
the following we will formulate our results for this transformed random evolution 
equation (^). 

Note that ri — G{hrj) can be handled as a function fulfilling the boundary conditions 
with A(77 - G(A?7)) 0. 

3. Existence and uniqueness of the solution for the 3D 

QUASIGEOSTROPHIC EQUATION 

In this section we consider the well-posedness of the transformed QGE (^. We 
intend to consider equation (|^) as an evolution equation on a rigged space V C 
H d V' . Since this equation is given with respect to periodic boundary conditions 
in horizontal directions the space H is defined to be 

L2{0) = {uEL2iO), [ udO = 0}. 

Jo 

The inner product on H is denoted by (•,•)■ Let H'^, A; G N be the usual Sobolev 
space consisting of functions with square integrable derivative up to A:-th order. If 
k is not an integer these spaces are the usual Slobodeckij spaces. For V C n H 
we choose the space of functions from such that 

Note that for functions in the trace on the boundary is well defined. This set 
will be equipped with the usual inner product denoted by (•, Let ^ be a 
linear bounded operator 

A:V 

which is the usual operator stemming from the bilinear form a{u, v) defined by —A 
with periodic boundary conditions in x,y directions and homogeneous Neumann 
boundary condition in z direction. Note that —A is symmetric with respect to 
these boundary conditions. We have that 

{Au,u) = 
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(•,•) denotes the dual pairing between y = D{A^i) and y = D{A^). The operator 
A is an isomorphism from V to V' . Since V is compactly embedded in H we can 
define G(/) = A~^f which is a continuous operator for instance from H to D{A) 
or more generally from Z)(A*) to D{A^^^), s G M. For the definition of the spaces 
D{A'^) and their norms || • ||_d(a=) see Temam |30) Section II. 2. In particular, for 
f G H the function G{f) is periodic in x, y directions. 

On H the operator A has the spectrum < Ai < A2 < • • • . The associated 

eigenf unctions are denoted by ei, 62, • • • . 

The embedding constant between V and H is given by Ai: 

^i\\ufH<\\u\\l. 

We now investigate the properties of the operator J which defines the nonlinearity 
of®. 

Lemma 3.1. Suppose that mi, m2, TO3 are three nonnegative numbers less than | 
with 

' 3 
2^m, > -. 

1=1 

Then there exists a constant ci > such that for u G H™'^^^ , v G iJ^s+i ^^^^ 

\{J{u,v),w)\ < ci||u||^™i+i||w||h™2+i||w||h".3. 
Proof. On account the Holder inequality, we get 

\{J{u,v),w)\ < ||V«|U,^ W^vUJuU,^, 1 + 1 + 1 < 1. 

' ' ^ qi 92 qs 

The embedding property iJ™ C Lq{0) where | = i- ^forg>l,m>0, to^| 

gives the conclusion. For the idea of the proof see Temam [^, Section 2.3. □ 

Remark 3.2. Suppose two of the m^'s in the last lemma, say m2, m^, have the 
value zero. Then if toi is chosen bigger than | the conclusion of the last lemma 
remains true. This follows if we apply at first the Sobolev lemma and then the 
Cauchy-Schwarz inequality. 



From Lemma 3.1 we can derive some a priori estimates for the nonlinearity of (Q). 



Corollary 3.3. Suppose that u G D{A°'), a > | and w . Then 
\{J{u,v),w)\ < ci\\u\\d(A'')\\v\\v\\w\\v , 

and hence for u G H, v, w Cz V 

\{J{G{u),v),w)\<c,\\u\\h\\v\\v\\w\\v. 

Suppose now that u G D{A°'), a > | and v e V, w e H. Then 

\{J{u,v),w)\ < Ci\\u\\D(A^f\\v\\v\\w\\H, 

and hence for u G D{A"), a > j and v E V, w G H 

\{J{G{u),v),w)\ < Ci||M||£,(^»)||u||y||w||//. 

Furthermore, if u E and v, w V we have the inequality 

\{J{U,V),W)\ < Ci||M||//2||w||y||w||v/. 
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Based on this Corollary, we also get 

Lemma 3.4. For any > there exists a constant C2,^ such that for v ^ H and 
\{Jiu,v),Av))\<c,J\u\\UMl + f^\\MH- 

and for an rj E 

\{J{v,v),Av))\<c',Jrj\\H4v\\^^^l^\\Av\\o< 

For the last inequality we need an interpolation argument and the Young inequality. 
Later on we can use this lemma to derive additional regularity properties of (Q). 
Also on account of Lemma 3.1 we obtain the following algebraic properties of J. 

Lemma 3.5. Ifv, w eV andu G and periodic inx, y directions, then we have 

{J{u,v),w) = -{J{u,w),v). 
Hence we have {J{u,v),v) — 0. 

Proof. We choose u from a set of sufficiently smooth functions which is dense in 
D{A"). Integration by parts formulae yields 



UxVywdxdydz — I I UyVxwdxdydz 



Uxyvwdxdydz + / Uyxvwdxdydz 



2-iT r /■It; 

UxVWydxdydz + / UyVWxdxdydz 

UxVw\ _p dxdz ~ / / UyVw\^_^ dyaz. 
^~ Jo Jo ^~ 

Note that the two boundary terms are zero. Indeed, since u(x,0,z) — 27r, z), 
we know that Ux{x, 0, z) = Uxix, 2tt, z), and since v, w are 27r-periodic with respect 
to X. By the smoothness assumption we can suppose that the derivatives on the 
boundary are well defined. Thus the first boundary term is zero. Similarly, the 
second boundary term is also zero. By the continuity of J (see Corollary 3.3) we 
can extend the relation to x V x V. The second claim is due to antisymmetric 
property of the first claim of this corollary. □ 

The following lemma will be used to obtain the continuity of the solution operator 
of®. 

Lemma 3.6. There exists a constant C3 > such that for ui, U2 E V 

\{J{G{ui),Ui) - J{G{U2),U2),UI - 1*2)1 < C3IIM1 - U2||y||ui - U2||H||ui||y. 



Proof. By Lemma 3.5 the left hand side above expression is equal to 

|(J(G(ui),1ti),M2) + (J(G(U2),U2),U1)| 

= |(J(G(ui),Ui),U2 - Ml) - (J(G(u2),Wl),M2 " Wl)| 
= \{J[G{ui) - G{u2),Ui),Ui - U2)\. 

Corollary ^.3| gives the conclusion. □ 
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The equation can be considered as an evolution equation on the rigged space 
V C H C V introduced at the beginning of this section. 

We now explain the properties of coefficients which are contained in (Q). Due to 
Corollary [3.3| we have a bilinear continuous operator 



B{-,-) -.V xV ^ H, B{u,v) := J{G{u),v). 

In addition, we have a time dependent linear operator 

C{t, u) J{ri{t) - G(A7y(t)), u), Cit) -.V^V 

where ri{t) G L2 ;oc(0, 00; HH) and periodic in x, y directions. Then G{Ari{t)) e 
DiA). 

Lemma 3.7. Suppose that r]{-) G L2,;oc(0, cxd; and 77 is periodic in x,y direc- 
tions. Then we have for almost any i > 0; 

\\C{t,-)\\vy. <CiMt)\\H2, {C{t,u),u)=Q. 

Indeed, because rj{t) — G{Ari{t)) G for any t > we obtain the first conclusion 



by Corollary tO. The second Conclusion follows from Lemma 3.5 



We now investigate the last linear operator appearing in which is defined by 

Diu) = (3G{u)^ -.V-^H. 
We obtain straightforwardly ||-D(-)||v.ff < C5. 
Lemma 3.8. For u E V we have 

{D{u),u) = 
Proof. Denoting G{u) G D{A) by ^ we have 

f3-HD{u),u) AO = ]^l^jj,C)xdxdydz - {ey)xdxdydz 

F{z) / iOxdxdydz 



Jo 

via the integration by parts. The second term under the integral generates the 
boundary term 

7: / / ^xiv\^ J' dxdz. 
2 Jo Jo 1 2^=0 

Indeed, for sufficiently smooth ^ from a dense set in V we have ^(x, 0, z) = ^(x, 2tt, z) 
and thus i,a:{x,0,z) = ^xix,2'K, z). It follows from the periodicity in y we have 
^y{x, 0, z) = $,y{x, 2tt, z) such that this boundary term is zero. 
For the last term the following boundary term appears 



2tt p2tt 



z=2tt 

dxdy 

z=0 



which is zero by the homogeneous Neumann boundary conditions. Integration with 
respect to x and using the periodicity as in the proof of Lemma 3.5, we get the 
conclusion. □ 
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We can formulate the evolution equation on y C i? C y 

(5) ut + i^Au + B{u,u) + C{t,u) + D{u) = f{t), u{0) = uq e H 

for f{t) — [3{G{Ki-j{t))x — rj{t)x) which is contained in L2,ioc{0, oo; H). 

Apart from the linear operators C {t) and D equation (||) has the form of equations 

of 2D Navier -Stokes type; see Temam pS], Chapter 3, for which we have existence 



and uniqueness. Here 2D means that the conclusion of Lemma is fulfilled 
which is responsible for a uniqueness theorem. By the regularity properties of the 
operators B, C{t) and D the same method as for the 2D Navier-Stokes equations 
ensures existence and uniqueness for (^. Thus we get the following main result in 
this section, about the well-posedness for 3D quasigeostrophic flows under random 
wind forcing on ocean surface. 

Theorem 3.9. (Well-posedness) Suppose that the Ornstein-Uhlenbeck process rj 
defined in is in L2,ioc(0, cxd; Then the 3D QGE equation or ^ has a 

unique solution u{t) £ -/j2,ioc(0, oo; V) D C([0, oo); H). 



Remark 3.10. i) On account of Lemma 3.6 we can prove that u{t) depends con- 



tinuously on the initial conditions x E H . This will be used later on. 



ii) Due to Lemma 3.4 for t > and a bounded set of initial conditions in H the 
image with respect to the solution operator uo u{t) is bounded in V . Hence the 
solution operator is compact for i > 0. 



4. The dynamical system of 3D quasigeostrophic flows 

In this section we study the dynamical behavior of QGE (||). In the following we are 
going to describe the background perturbations defined on the boundary of O which 
will influence the dynamical system generated by (^. We will have two different 
influences. The first perturbation is a white noise which models the weather or the 
small scale impact of the atmospheric motion through wind forcing on the surface. 
The other one is a periodic motion which serves as a model for the impact due to 
periodic rotation of the earth and thus periodic exposure of the earth to the solar 
radiation; see Chapter 6 and Chapter 11. 

In the first part of this section we are going to explain a dynamical model of the 
boundary conditions. 

We consider the elliptic differential equation 

At. = 0, 1^ = on a...o 1^ = / e H^+"^ {0,,2.) n L2(0....2.). 
on on 

On the other faces we have periodic boundary conditions. This equation has a 

unique solution; see Egorov and Shubin Q Page 130f. The solution operator is a 

linear continuous operator with the image in H^'^'^ . We denote this operator by 

G(/). (O. . 27r) denotes a usual boundary space. 

At first we consider the random part of the boundary conditions. Let be a con- 
tinuous temporal Wiener process with values in a Hilbert space U — i?^'^^ (O. . 2^)n 
-^2(0. ^.^27r) This Wiener process is defined for positive and negative times; see 
Arnold ^ Page 547. The covariance operator of W is denoted by Q, which is 
a positive symmetric linear operator on U with finite trace. The dynamics of W 
is given by the metric dynamical system consisting of a probability space (fi, P) 
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and a flow 9, (ri,jr, P, 6*) where P is the Wiener measure with covariance Q and 
^ = (^'t)teR is the flow of the Wiener shift: 

W{-,9tUj) = W{- + t,uj) -W{t,uj) foriGR. 
The mapping 9 is (B(R) ® T ^ JF) -measurable and fulflUs the property 

(6) 9t+r=Gt°Qr, i, reR. 

For instance, we can choose ^2 to be the set of continuous functions Co(R, U) which 
are zero at zero and for T we choose the Borel-cr-algebra of Co(R, i7). Note that 
P is ergodic with respect to 9. Later on we have to restrict this metric dynamical 
system to fix particular dynamical properties. 

We now show the existence of a solution of (^) satisfying particular properties. At 
first we show that the following problem has a stationary solution: 

(7) mt^^m, ?r^-0, e O.,.,o, ^^W{t), {x,y)E0,,2n 

on on 

under periodic boundary conditions with respect to the other faces of the cube O 
and a initial condition at time t = 0. This solution will serve as a process which 
compensates the nonhomogeneous boundary conditions in (^. A similar problem 
has been considered in Da Prato and Zabczyk |^ Chapter 13 or p^ . 
In contrast to (|^) we now consider boundary conditions which are defined to be ho- 
mogeneous Neumann boundary conditions on 0.,.,27r- The solution of ^ with these 
boundary conditions generates an analytic Co-semigroup (S(t))t>o. This semigroup 
has a generator, which will be denoted also by —vA where A is equivalent to the 
operator introduced in the last section. 

The following lemma allows us to define a stationary Ornstein-Uhlenbeck process 
which fulfills particular regularity assumptions. We notice that G maps continu- 
ously Hi{0{-, •, 27r)) into V H where 

V^iueH^nH : u(0, •, •) = u(27r, •, •), 0, •) = 27r, •)}• 

We also note that for our boundary conditions V fl ^ D{A). 

Lemma 4.1. Let Q be a linear operator of finite trace onU ~ Hi (0.,.,2-n)^H ■ We 
consider the elliptic differential operators with respect to the boundary conditions 
formulated in 

Av + 2{F^v)^ + F,,v + F,,,u 

which are supposed to be generators of analytic Co-semigroups of negative type. 
Then we have 

\\AS{T)\\l^(u,m^ < oo, / < oo, /3 < -. 



Proof. Let be a complete orthonormal base in U . Then GQ^Ci G H^C]H. We now 
study the D^, D^x, ■ ■ ■ i^zz-derivative of GQ^ei. All these expressions are periodic 
in X, y directions. In addition, all these elements are contained in D{A^), ?7 < | 
This follows because we only have periodic and Neumann boundary conditions, see 
Da Prato and Zabczyk |^ Page 401. Since S{t)g is the solution of a differential 
linear parabolic equation 

(9) ^ + ^Av = 0, v{0) = g. 



(8) Au+{F,u),, 
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we can differentiate the the solution with respect to D^- The coefficients of the 
equation are independent of x. The function DxS{t)g fulfiUs the same equation as 
S{t)g but with derived initial condition. Hence we have DxS{t)g = S{t)Dxg. The 
regularity of D^g G D{A^) ensures that \\ADxS{t)g\\H is well defined such that 
D^S(t)g e D{A), see |l| Page 392: 

ll^xyi\D{Av) 



(10) \\AD,S{T)g\\l^\\ASiT)D^g\\l<-^e~'"'\\D,g\\l.^^ t > 0. 



for g = GQ^ii, a > and rj £ (|,|)- I3ut on D{A) we can identify A with 
—A such that we can exchange A and D^- Integration of (10) from to oo and 
summation over i yields the desired estimate. Similarly, we can argue for the 
derivatives Dx , Dy , Dyy , • • ■ . 

Similar considerations allow us to prove the estimates for derivatives Dzx, • • ■ , Dzz- 
In contrast to the above formulae a derivative Dz changes A to an operator A given 
in (||) with semigroup S{t). The assumption (^ ensures that S{t) is of negative type 
with our usual boundary conditions. We get 



\\DzAS{T)gfH = \\ADzS{T)g\\l = \\AS{T)DzgfH < -^^-^^'\\D.g\\\j,,y t > 0. 

The same method gives us estimates in z-direction. Note that if we differentiate 
the coefficients of (^ twice in z direction one term without a second z derivative 
appears. To take this term into account we need the system of differential operators 
in (I). 

To obtain the second conclusion we use the estimate for /? e (0, |) 

\\ASit)g\\l^^,^ = \\A^ASit)g\\l < -^e-^'^'MliA^^^.y 

and an appropriate 77 > i; 77 + /3 e (i, |), t > 0. □ 

Note that if N{z) (hence F{z)) is constant then the assumption about the generator 
is superflow. 

Remark 4.2. On account of Da Prato and Zabczyk |^ Chapter 13 there exists 
an Ornstein-Uhlcnbcck process 

m{;^)^A f S{--T)dGW{T,uj)€L2,ioc{0,'^;H^) 
Jo 

solving (|^) with initial condition zero at t = 0. ?7i(i, ui) is periodic in x, y directions. 

Now we are going to show that r]i{-,Lu) has continuous paths in H. 

Lemma 4.3. There exists a continuous version in H and in particular, 

E sup ||?7i(s)|||^ < 00. 
se[o,i] 

Proof. Due to Da Prato and Zabczyk [|l] Theorem 5.9 and Remark 5.11 we have 
to check 



POO 

/ T-^||A^(T)GQi||2 dr<oo 
Jo 
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for some 7 > 0. The left hand side is equal to 
(11) 



see Da Prato and Zabczyk ||2^ Proposition 13.2.4 where (e^) is orthonormal base 
introduced above. For sufficiently small 7 > we have a 0(7) > 0: 

Jo Jo 
Hence ( |lT| ) can be estimated by 
00 „t 

00 I't pt 

<c{j)J2x^+'^\\QiG*e,\\j, e-2A.-dr = c(7) / \\AS{r)\\l^^,,^^,^^dr 
i=i "^0 -^0 



which is finite by Lemma 4.1. □ 



We now are looking for a stationary process solving (|7|). 

Lemma 4.4. There exists a random variable rji G H such that t — s- rjiiOtUj) solves 
#• 

Proof. By the Fourier method we have ||<5'(i)|||/^ < e"^'*'^*. We obtain for ni < 
77.2, ni ^ 00 

E sup ||77i(s + Til, 6'_„iCj) -771(5 + ^2, 6'_„2Cj)|||j- 

s6[0,l] 

rni+s 

(12) =E snp \\A S{ni+s-T)dGW{T,e^n,uj) 

^ ' se[oa] Jo 

I 2 



-A S{n2 + s-T)dGWiT,0^n,co)\\j, 

Jo 

We can continue by a simple integral substitution, the ((?4)fgR-invariance of P and 
the semigroup property of 5 

= E sup \\A f S{s~T)dGW{T,uj)~ A [ S{s~T)dGW{T,Lo)\\]j 

sG[0,l] J-rii J -112 

= E sup \\A f ' S{s-T)dGW\\% sup \\A f S{ni + s - T)dGW\\jj 

sG[0,l] J-n2 se[0,l] J~n2+ni 

< sup e-2-^i("i+^)E||A / S{^T)dGW{T,Lu)\\jj 

s€[0,l] J -ni+n2 

2 ^ ^ — 2Aini 



/>oc 

<e-2^^"^ / II A5(r)GQi ||2 <C7e 
Jo 



The finiteness of the last integral follows by Lemma O Hence 771 (■ + 77, 0_„w) is 
L2-converge nt w ith respect to C([0, l];-ff). This convergence is exponentially fast. 
Via Lemma 4.3 and the semigroup property the exponential convergence to zero of 

S{s-T)dGW{T,^)f^ 

-11 
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for n — > cxD is easily seen. 

By the Borel and Cantelli Lemma there exists a set fi^ C of full measure such 
that for w € we have the above convergences w-wise. The same remains true on 
fJ-' if we replace w by OjUi for j G Z. The set fl^ is (6't)tgz-invariant and has a 
full measure. 

We denote by rji {uj) the w-wise limit of rji (0 + n, 0_„w) which exists on . 
To see that 771(0 + n, is converging on we can restrict ourselves to 

t G [0, 1]. In particular, if f G (— 1, 0), B^uj can be written as 9i_^_t6^iLLj, 1 + f G (0, 1) 
for i G (—1, 0). We can write 

771 {n, e^^+iuj) = A f S{t - T)dGW{T, oj) 

/i p — n~\-i 

S{t-T)dGW{T,Oj)- A / S{t~T)dGW{T,Oj). 
- n J —n 

By the above considerations the first term on the right hand side is convergent 
and the second term tends to zero by ( p^ using the semigroup property. Let Cl 
be the set of oj's such that rii{n,9-nUj) converges w-wise. We have for i G M that 
6*4^2 D pl^ il^ . Hence the (6'()(gK-mvariant set HteR^*^ contains a set of measure 
one. We can construct a new probability space if we consider the trace of J- and 
the restriction of P to this new set f]^ 6t^. We will use for this new probability 
space the old notations. 

To see that t rii{9tuj) solves we consider rii[6tLo) — S{t)rii{uj) which is a version 
of A /J S{t — T)dGW such that by the variation of constants formulae 

s{t)'ni{uj) + A I s{t-T)dGW{T,uj)^'qi{etoj), 



see Da Prato and Zabcyk [|| (13.2.13). □ 

Remark 4.5. Using the technique for the proof of the last lemma we can also show 
that there exists a (^^t)^gR-invariant set of full measure such that t — » rji^OtUj) has 
paths in L2,ioc(0, c»; H^). Since the distribution of rji is a Gaufi measure it can be 
derived that the stationary process 771 (^.) has trajectories in i4,ioc(0, 00; H"^ Ci H), 
see Da Prato and Zabczyk Corollary 2.17. 

We now introduce the term tempered random variable. A random variable ^ is 
called tempered if this random variable has a subexponential growth: 

y log+(C(g«a;)) „ 
limsup j-j — U. 

t-toc\teT \t\ 

In the case of an ergodic measure P, it is well known that there is only one alterna- 
tive, i.e., the limsup in the last formula is +00; see Arnold ^ Page 165, Proposition 
4.1.3. 

Lemma 4.6. There exists a {Ot)t^s.- invariant set VLi of full P-measure such that 
for any lo ^ Qi . 771 (cj) is well defined and 

l0g+||77l(g,^)||ff _^ 
t^±oo t 

which means that the H norm of 771 is tempered. 
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Proof. On account of (12) we have 

E sup \\rii{9sUj)\\H < oo. 

sG[0,l] 

Hence by the Birkhoff ergodic theorem, we have a (0f )tgz~invariant set fJi C 
such that 

SUp,g[o,l] \\Vl{Ss+n^)\\H ^ 

hm = 0. 

On the other hand for t > 0: 



such that 



\m{^tUj)\\H < sup \\rii{es+[t]Uj)\\H 
se[o,i] 



hm = 0. 

t^±oc t 



and similarly for t < 0, which is sufhcient for the (0t)igR-invariance of fii. □ 
Collecting our calculations we get 

Theorem 4.7. Let Q he a linear operator on U — {0.,.,2tt) H L2{0.,.,2tt) with 
finite trace. Then there exists a tempered random variable rji and a (9t)teR-'i'™o,i^io,i^t 
set of full measure such that 

t -> 7ni9ti^) e C([0, oo); ff) n Li^odO, oo; n H). 

r]i{Lu) is periodic in x, y directions. 

In a similar manner we can consider (^ with time-periodic boundary condition 
representing the impact of the earth's rotation on the fluid : 

(14) ?72t = ^V2, -T^ = 0, on O.^.^o, = wo sin(27rt) on 0.,.,27r 

for Mo G n L2{0....2tt) and (spatial) periodic boundary conditions on the other 
faces. We obtain without proof the following : 

Lemma 4.8. Suppose that uq G H^{0. ...2Tr) n i2(0.,.,27r). Then there exists a 
continuous periodic solution 

t^ri2{t)eH^ 

which satisfies ^I^. In particular, 772 (i) also periodic in x-y directions. 
The proof is based on the properties of the operator G. 

Let 9"^ — {9'^)teR be the shift operator 9lf{-) = f{t + •) for appropriate functions 
/. We consider the hull of uq sint with respect to 9^: 

^2 = U ^'?(wo sin(27r-)) ^ [j uq sin(27r(i + ■)) = |J sin(27r(t + •))■ 

teR teK te[o,27r) 

Summarizing, we have found a process ry = ?7i + 772 which will serve as a model for 
the perturbation on the ocean surface. 



After these preparations we can introduce a nonautonomous/random dynamical 
system. Let 6* be a flow on a set fl (such that (||) is fulfilled). 
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Definition 4.9. A nonautonomous dynamical system (f) on a phase space H with 
respect to a flow 6 is a mapping 

0: T+ xflxH ^ H 

fulfilling the cocycle property 

(j){t + T, OJ, ■) = (t){t, OrUJ, (j>{T, UJ, •)) for t, T e T+ 

0(0, t^, •) = iAh 

for a; G and t, t G T+. Suppose that the flow 9 is carried by a metric dynamical 
system and is supposed to be measurable then (p forms a random dynamical system. 



This definition is due to Arnold or Kloeden et. al. |15 



Remark 4.10. We now consider the flow 6 — {9^,6'^) on = fi^ x fi^. By 
the global forward existence and uniqueness of the solution of (H) for ri{9tuj) — 
r\\(9\L0\) + 7/2(0^^2), ^ = (wi, W2) G f^, the solution operator of (^), which maps an 
initial condition uq G H and a sample point lo to the solution at time t, has the 
cocycle property for H = L2{0). We will denote this solution operator by x). 
Note that by the periodicity of 772 the restriction of on Z is a random dynamical 
system for any fixed uj2 G ^2- Indeed, 9iUj = {9luji, UJ2) which can be identified with 
9^ and which leaves P invariant (P is ergodic with respect to {9})i^z,)- 



Remark 4.11. Another opportunity to get an example for a complete random 
dynamical system would be to equip f22 with an ergodic measure. But in contrast 
to the fact that the daily or yearly rotation of the earth is well determined such a 
random ansatz would express that the beginning of these periods is rather random. 

The main result of this article is the existence of an attractor for the nonautonomous 
dynamical system. This attractor will attract random sets in probability. Before 
we give the main theorem, we make some basic remarks on random sets. 
Suppose that H is a. Polish space. A set function u; — *■ D{uj) with closed and 
nonempty images is called a closed random set over (f2i,jr, P) if and only if there 
exists a countable number of random variables 

: (r!i,^,P) ^-ff, leN 

such that 

see Castaing and Valadier Q Chapter 3. 

A random set D is called tempered if the random variable dist//(_D(ijj), {0}) is 
tempered where 

dist^l^, B) ~ sup inf \\a — b\\H. 

aeAbeB 

We now define the term random attractor. 

Definition 4.12. Let <j) be a random dynamical system over the metric dynamical 
system 9. A tempered random set (yl(aj))t^gsi with compact and nonempty images 
A{lo) is called random attractor if 

(15) (j){t,w,A{uj))= A{9tLo), <eT+ 
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and for any random tempered closed random set D{lj) 



(16) (P) lim distH(Mt,uj,D{uj)),A{etuj)) = Q. 

In the next section we will show that the dynamical system generated by re- 
stricted to T = Z has such a random attractor. However, because of time-periodic 
perturbation, we do not have a random dynamical system but a nonautonomous 
dynamical system for T = M. Therefore we have to modify this conclusion for 
T = Z a little bit, and thus we obtain the following main result in this section. 

Theorem 4.13. There exists set function A = {A{uj))^^q, for the 3D QGE M) 
under random plus time-periodic forcing on the ocean surface such that ( |7^ j (|7qjzg 
satisfied for set functions D with closed images and subsexponential growth of t ^ 
dist//(I?(0tw), {0}) with tempered restriction to {9t)tez- Moreover, the mapping 

t ^ AiOtLj) 

has a subexponential growth. 

In the next section, we prove this main theorem 4.13| . 



5. Proof of the main theorem 



The proof of the main theorem 4.15 in the last section is based on checking following 
conditions. We will consider the set of families of sets D such that D{u!) is closed 
and nonempty and t ~> dist h {D (0 t{uJi,LU 2)), 0) has a subexponintially growth. If 
we consider D with respect to the restriction of {9t)t£M to Z family of sets D is 
supposed to be tempered. The set of these set families is denoted by V. 

The following theorem states the existence of random attractors. 

Theorem 5.1. Suppose that 6^ be a metric dynamical system and suppose that for 
LVi G ril, t e Z+ the mappings 4>{t,UJi,-) are continuous and that there exists a 
tempered set B having compact images which is absorbing: 

(17) 0(t,0Vi,^(^Vi)) cB(wi) 

for Z+ 3 t > to{D,uj), D is supposed to be tempered. Then there exists a unique 
random attractor A which is tempered. 

The proof of this theorem can be found in Flandoli and SchmalfuB [0 . 



Now let us use this result to prove our main Theorem 4.13 



Proof. Calculating the inner product in H we obtain by Lemma 3.5| , 3.7, 3.i 
(18) 

\\u{t)fH + 2iy f \\u{T)\\UT=\\uo\\l + 2f3 f {G{Kr,), - r,,,u)dT 



<\\uo\\% + — I \\G{Ar,{erU;)U-v{erU;),\\l,dT + iy I \\u{T)\\l.dr. 
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Parallel to this inequality we consider the affine equation 

-||«0|Ih + — / \\G{^T^{erU0))^~T^{erUjU\l'dT. 



The solution of this equation ^{t, uj, ||wo|||/) is a bound for |||^. In addition, this 
equation has a unique forward and backward exponentially fast attracting solution 

{cu,t) -^CiOtio), rH = — r e''^'^\\G{Af^{drUj)h^v{OrCoU\v'dr. 



for t —>■ oo which follows by the variation of constants formulae. 

To see the backward convergence uniform with respect to D we must prove that 

hm sup |C(t,0-ta;,ho||l/)-rMI =0. 

Because of 
we have 

sup |C(t,0_tc^,ho||lf)-rHI < sup |||^/o|||^-r(^-t^)|e-^^* 

•ito6-D(f_ti^) uo6-D(e-ti^) 

which also follows by the variation of constants formulae. Since the first factor on 
the right hand side is only subexponentially growing the convergence conclusion 
follows for f — !■ cxD. 

The forward convergence follows straightforwardly by the variation of constants 
formulae 

I sup at,^AWofH)~CiOtto)\<e~^'' sup \\\uo\\l-eH\ 

where £^*(9tLu) = S^{t,LU,^*(Lu)). The mapping t — > S,*(6tuj) is subexponentially 
growing. Indeed, we have 

cO 



— oo 




The second term on the right hand side is bounded and the first is tempered, see 
Theorem 4.1 in H^. Since ||77i||h is tempered so is ||G'(A77i)2;|| v 
Hence the ball B with center zero and radius 2^*{uj) forms an absorbing set (in the 
sense of (|l|)), B{lo) G V. Plugging in the radius 2^*{lu) of 13 into ® for \\uo\\jj 
it is easily seen that B is forward invariant: 

<j>{t,uj,B{Lu)) C B{9tLj). 

Since for fixed lu2 the mapping ^* is a random variable B is a random set. On the 
other hand we can check that 



B{lu) := <l){l,9^iUJ,B{e^iLu)) C B{uj). 



Note that B is compact because (/> is regularizing which follows by Remark 3.10 



ii). We have seen that B is forward and backward absorbing. This remains true 
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if we restrict our random dynamical system to discrete time, see Remark 4.10 
By the above inclusion we can derive B G V since B G V. It follows also by the 
definition of B that for fixed [02 the set B{(uji, L02)) is random set which is tempered 



The continuity of w, •) follows because we have the estimate from Lemma 3.6 



However, we need this technique later on once more such that we are going to 



explain this technique: We have by Lemma 3.6 and the properties of C{t) and D 



\{B{ui,ui) - B{u2,U2),ui - 7/2)1 < C3IIM1 - U2\\h\\ui - U2||yl|uil|v 

{C{t,Ui)^C{t,U2),Ui~U2) =0 
{D{ui) - D{U2),UI-U2) -0 

Hence there exists a constant cg such that for two solutions ui, U2 with initial 
conditions uio, U20 

- H2WIII, < C6(||t.iWfy)||t^lW - ^2(t)|||^. 
It follows by the Gronwall lemma 

(20) \\ui{t)-U2{t)fH < \\uiO-U2o\\le^'^o\\MrnUr 

which gives the continuity of (j){t, lu, •). 

We have checked all assumptions of the above theorem such that the discrete ran- 
dom dynamical system has a random attractor A{uji,uj2) for some fixed lu2 with 
respect to the restricted discrete random dynamical system. We now extend the 
definition of A to n2. We set: 

A{0t{ui,Q2)) = Ht, K,(I'2),v4(wi,W2)) for t e K+ 

which is invariant in the sense of (|l^). Note that by the cocycle property this 
definition is correct. We show convergence (fiW: 



(P) lim distff(0(t,(u;i,(52),£'((wi,W2))),A(0tK,W2))) =0. 

for D €z T) and any cj2 G ^^2- Since B(uj) is a forward absorbing and forward 
invariant set for any _D G P it remains to check the convergence conclusion for 
D = B. On account of (|l|) we can notice for fixed uj2 and ui ~ (j){t, (wi, (^2), y) 

sup V I ||0(t, (a;i,aj2),y)llvrfr 

yG-B((wi,i32)) "'0 

<2rK,cS2) + — / \\G{^'n{er{uJi,UJ2))y^-fl{6r{uJl,id2)yA\l'dT. 

With respect to we have an appropriate constant c^ > 0: 

sup e^e Jq ll0(T,(a.i,a>2),a:)|lv<i'^ 
x^B{{u)i,a2)) 

< gCr(5*(c^i,i'2)+Jo l|G(Ar,(e,(a;i,i>2)):„-»7(ex((^i,<S2))x||^,dT) _ Y{{uJi,Lj2)) 

for t G [0, 1]. The mapping 



THE STOCHASTIC 3D QUASIGEOSTROPHIC FLOWS 



19 



defines a stationary process. We obtain by (^0|): 

sup dist/f-(0(T + n, UJ, (Ii2)), A(e'r+n(t^i, CJ2))) 

re[o,i] 

< distif(0(n,(wi,(I>2),B((wi,<I>2)),A((6'iwi,(I'2)))l^((6'>l,W2))i 

Since the first factor of the right hand side tends to zero in probabihty the product 
of the right hand side also tends to zero in probabihty which gives the convergence 
conchision (p^. 

To obtain the subexponential growth of 

Ai0t{LUi,LU2)) 

we need that Y is tempered with respect to the restricted fiow {dt)tei,- Due to 
Arnold Proposition 4.1.3 it is sufficient to show that 

E^*(wi,W2) < 00, E \\G{Ar]{er+s{cJi,u'2)))x - r]{(^T+s{oJi,uJ2))x\\v'dr < 00 
Jo 

which can be derived from Lemma 4.3. Consequently, 



sup \\x\\h <Y{{LUi,LU2))''diam{A{{LUi,uj2))) + sup \\x\\h, te[0,l] 

by the triangle inequality. We have used that the product of two tempered random 
variables is tempered which gives the general convergence conclusion of the main 
theorem 4.13| since A{lui,uj2) is (6't)tgz-tempered. □ 



6. Summary 

We have studied the 3D baroclinic quasigeostrophic flow model under random wind 
forcing and time-periodic fluctuations on fluid boundary, i.e., on the interface be- 
tween the oceans and the atmosphere. The time-periodic fluctuations are due to 
periodic rotation of the earth and thus periodic exposure of the earth to the solar ra- 
diation. We have established the well-posedness of the baroclinic quasigeostrophic 



flow model in the state space (Theorem 3.9), and we have demonstrated the exis- 



tence of the random attractors (Theorem 4.13), again in the state space. We have 



also discussed the relevance of our results to climate modeling. 

Acknowledgement. A part of this work was done at the Oberwolfach Math- 
ematical Research Institute supported by Volkswagen Stiftung while the authors 
were Research in Pairs Fellows. This work was partly supported by the NSF Grant 
DMS-9973204. 



References 

[1] L. Arnold. Random Dynamical Systems. Springer, New York, 1998. 

[2] L. Arnold. Hasselmann's programm visited: The analysis of stochasticity in deterministic 

climate models. Report 450, Universitat Bremen, Institut fiir Dynamische Systeme, 2000. 
[3] A. J. Bourgeois and J. T. Beale. Validity of the quasigeostrophic model for large-scale flow 

in the atmosphere and ocean. SIAM J. Math. Anal, 25:1023-1068, 1994. 
[4] J. R. Brannan, J. Duan, and T. Wanner. Dissipative quasigeostrophic dynamics under random 

forcing. J. Math. Anal. Appl, 228:221-233, 1998. 
[5] C. Castaing and M. Valadier. Convex Analysis and Measurable Multif unctions. LNM 580. 

Springer- Verlag, Berlin-Heidelberg-New York, 1977. 
[6] B. Desjardins and E. Grenier. Derivation of quasigeostrophic potential vorticity equations. 

Adv. Diff. Eqns. 3 (1998), 715-752. To appear. 



20 



JINQIAO DUAN AND BJORN SCHMALFUSS 



[7] The dissipative quasigeostrophic equations. The dissipative quasigeostrophic equations. Math- 

ematika, 18:265-288, 1981. 
[8] J. Duan. A remark on three dimensional baroclinic quasigeostrophic dynamics. Appl. Math. 

and Comput., 106:285-289, 1999. 
[9] Yu. V. Egorov and M.A.Shubin. Partial Differential Equations, volume 1 of Encyclopaedia of 

Mathematical Sciences. Springer, New York, 1991. 
[10] P. F. Embid and A. J. Majda. Averaging over fast gravity waves for geophysical flows with 

arbitrary potential vorticity. Comm. PDEs, 21(619-658):619-658, 1996. 
[11] F. Flandoli and B. Schmalfufi. Random attractors for the stochastic 3-D Navier-Stokes equa- 
tion with multiplicative white noise. Stochatics and Stochastics Reports, 59:21-45, 1996. 
[12] K. Hassclmann. Stochastic climate models, Part I. Tellus, 28:473-485, 1976. 
[13] G. HoUoway. Ocean circulation: Flow in probability under statistical dynamical forcing, in 

Stochastic Mod,els in Geosystems, S. Molchanov and W. Woyczynski (cds.). Springer, 1996. 
[14] D. D. Holm. Hamiltoiiian formulation of the baroclinic quasigeostrophic fluid equations. Phys. 

Fluids. 29:7 8, 1986. 

[15] P. E. Klocdcn, H. Keller, and B. SchmalfuB. Towards a theory of immerical dynamics. In 
H. Craucl and M. Gundlach, editors. Stochastic Dynamics. Springer Verlag, 1999. 

[16] M. Leroux. Dynamic Analysis of Weather and Climate. John Wiley & Sons, 1998. 

[17] P. MuUer. Stochastic forcing of quasi-geostrophic eddies. In P. MUUer R. J. Adler and B. Ro- 
zovskii, editors. Stochastic Modelling in Physical Oceanography, pages 381-396. Birkhauser, 
Basel, 1996. 

[18] C. Nicolis and G. Nicolis, Is there a climate attractor? Nature 311, October 11, 1984. 

[19] J. Pedlosky. Geophysical Fluid Dynamics. Springer Verlag, New- York, Berlin, 1987. 

[20] J. P. Peixoto and A. H. Oort. Physics of Climate. Springer, New York, 1992. 

[21] G. Da Prato and J. Zabczyk. Stochastic Equations in Infinite Dimension. University Press, 
Cambridge, Cambridge, 1992. 

[22] G. Da Prato and J. Zabczyk. Evolution equations with white noise boundary conditions. 
Stochastics and Stochastics Reports, 42:167-182, 1993. 

[23] G. Da Prato and J. Zabczyk. Ergodicity for Infinite Dimensional Systems. Cambridge Uni- 
versity Press, Cambridge, 1996. 

[24] R. Salmon. Lectures on Geophysical Fluid Dynamics. Oxford Univ. Press, Oxford, 1998. 

[25] R. M. Samelson. Stochastically forced current fluctuations in vertical shear and over topog- 
raphy, ,/. Geophys. Res. 94 (1989) 8207-8215. 

[26] B. Saltzman. Climatic systems analysis. Adv. in Geophysics, 25:173—233, 1983. 

[27] B. SchmalfuB. A random fixed point theorem and the random graph transformation. J. Math. 
Anal. Appl., 225(1):91-113, 1998. 

[28] R. Temam. Navier-Stokes Equation-Theory and Numerical Analysis. North-Holland, Ams- 
terdam, 1979. 

[29] R. Temam. Navier-Stokes equation and Nonlinear Functional Analysis. CBMS-NSF Re- 
gional Conference Series in Applied Mathematics. SIAM, Philadelphia, 1983. 

[30] R. Temam. Infinite-Dimensional Dynamical Systems in Mechanics and Physics. Springer- 
Verlag, Bcrlin-Hcidclbcrg-Ncw York, second edition, 1997. 

[31] S. Wang. Attractors for the 3D baroclinic quasi-geostrophic equations of large-scale atmo- 
sphere. J. Math. Anal. Appl, 165:266-282, 1992. 

(Jinqiao Duan) CORRESPONDING author: Department of Applied Mathematics, Illinois In- 
stitute OF Technology, Chicago, IL 60616, USA., Tel: (312)567-5335 Fax: (312)567-3135 
E-mail address, Jinqiao Duan: duan@lit.edu 

(Bjorn SchmalfuB) Department of Applied Sciences, University of Technology and Ap- 
plied Sciences, Geusaer Strasse, D-06217 Merseburg, Germany, 
E-mail address, Bjorn Schmalfufi: schmalfussSin.fh-merseburg.de 



